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The dynamics of many natural systems is dominated by non-linear waves propagating through the 
medium. We show that the dynamics of non-linear wave fronts with positive surface tension can be 
formulated as a gradient system. The variational potential is simply given by a linear combination of 
the occupied volume and surface area of the wave front, and changes monotonically in time. Finally, 
we demonstrate that vortex hlaments can be written as a gradient system only if their binormal 
velocity component vanishes, which occurs in chemical systems with equal diffusion of reactants. 


Introduction. Propagating non-linear waves are found 
in many extended systems, which can be of oscillating, 
bistable or excitable nature. Chemical examples include 
combustion [T], catalytic oxidation [2], and spiral pat¬ 
terns in oscillating chemical reactions [3]. In physics, 
dynamical phase transition boundaries are also known as 
domain walls [4]. In biological context, non-linear waves 
often fulfil an intra- or extracellular signalling role. They 
have been observed both on the micrometer scale, e.g. as 
calcium waves [5], and on a macroscopic scale for action 
propagation in neural [6] and cardiac [7] tissue, where 
they govern cognitive processes and coordinated contrac¬ 
tion of the heart. 

All above-mentioned systems share the emergent phe¬ 
nomenon of non-linear wave propagation. In contrast 
to classical electromagnetism, non-linear waves have a 
fixed amplitude and do not linearly superimpose when 
colliding; instead they may annihilate or reflect [8]. As 
a consequence, most non-linear waves (except solitons) 
lack the classical invariants of motion such as energy and 
momentum, which conservative systems do possess. This 
lack is no contradiction with the laws of thermodynam¬ 
ics, since generally energy is consumed in the process of 
wave propagation in the form of heat, which is why these 
waves are also referred to as dissipative structures [9]. 
Due to the lack of invariants of motion, the analytical 
and physical description of dissipative systems remains a 
challenge. 

Some progress has been made in rewriting the partial 
differential equations (PDFs) for specific dissipative sys¬ 
tems as a so-called gradient system m- Here, we are 
interested to write a general system as a continuous gra¬ 
dient system, which is a set of PDFs that can be written 
in the form m- 


±{s,t) = - — , (1) 


for a functional F = f{'K,ds^,s)ds. Here X is a state 
vector and s represents a set of parameters labelling a 
A^—dimensional spatial domain ft. The dot denotes the 
partial derivative with respect to time and the notation 
S./S. refers to the functional derivative [12]. 


Note that 0 immediately implies that 
dF f 5F ■ , f /SfY , 

(sx) ® 

Hence the potential F decreases in time and therefore is a 
Lyapunov function of the system. Thus, if a critical point 
is an isolated minimum of F, it will be an asymptotically 
stable equilibrium point [10]. In a sense, the dynam¬ 
ics of gradient systems is relatively easy to understand, 
since finding the potential function F at once reveals the 
long-term attractors of the dynamical system. However, 
finding the potential for a given system remains much 
of an art, and it is unlikely that a single method will 
be successful for any PDF that describes a dissipative 
system. 

We illustrate the difficulty of finding a gradient sys¬ 
tem form on the important class of reaction-diffusion sys¬ 
tems. These dissipative systems are modelled by a set of 
parabolic PDFs, which are in one spatial dimension of 
the form: 

u(x, t) = P9^u(x, t) + R(u(x, t)). (3) 

A necessary condition under which a system can be recast 
as a gradient system is found by taking the functional 
derivative of both sides in Fq. 0 . Letting the state 
variables act as the components of X in Q, we find 

SF Sii^ 

Su^ Su'^ 

This compatibility condition is the functional analogue 
of the well-known theorem in vector calculus that X = 
VF X X = 0 in a simply connected domain. For 

a reaction-diffusion system, the compatibility condition 
0 requires symmetry of the Jacobian of the reaction 
term. Thus, only for selected reaction terms R(u), one 
may rewrite (1^ as a gradient system, e.g. for the diffu¬ 
sion equation (R = 0) or when there is only one state 
variable m- With small amendments to the gradient 
system structure, the FitzHugh-Nagumo equations have 
been cast to this form [13] and extensions with skew- 
symmetry have been considered in M- However, con¬ 
dition Q implies that trying to rewrite every reaction- 
diffusion system as a pure gradient system is an impos¬ 
sible quest. 
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In this Letter, we take a different approach by for¬ 
mulating not the underlying system but the emergent 
wave front dynamics as a gradient system. Therefore, 
our method will work for any dissipative system that sup¬ 
ports traveling waves with uniquely selected wave speed 
and whose wave fronts end orthogonally to the medium 
boundaries. 

The dynamical law which we will formulate as a gra¬ 
dient system is the velocity-curvature relation for non¬ 
linear waves. In homogeneous media, the wave speed c, 
measured along the wave front normal n, is found to de¬ 
pend on the wave front’s extrinsic curvature K. When 
the wave front thickness is small compared to the front’s 
local radius of curvature, the velocity-curvature relation 
can be linearized to 


C = Co - (5) 

with plane wave speed cq and surface tension coefficient 
7 . For reaction-diffusion systems, the relation ^ ca n be 
derived from the underlying system equations [151 fT6| . 
By convention, the speed Cq is positive; it may vanish 
in certain dissipative structures m- Further, stability 
analysis shows that only when the surface tension 7 is 
strictly positive, the non-linear waves are stably propa¬ 
gating and we will work in this regime. In this regime, 
small protrusions of the wave front will be dampened, 
such that after a short while, the front’s curvature will 
be small, and it can be well described by Eq. (§. 

Theorem. Our main result in this work is that the 
geometric law of motion <§ for non-linear waves can be 
written as the continuous gradient system 0 with ever- 
decreasing potential 

F{t) = -coV{t)+^S{t). (6) 

where V{t) is the volume where the traveling wave has 
passed and S{t) is the wave front’s total surface area. 
Our result holds for wave fronts in any number of spatial 
dimensions N > 2. In what follows, we will use termi¬ 
nology of surface and volume as if working in N = 3 
dimensions. 

Proof. Our proof involves the calculus of variations 
and goes as follows. To describe the emergent dynamics, 
we first formalize the definition of a wave front. In a 
spatial domain that locally looks like we identify 
a scalar quantity u of the system (e.g. temperature or 
transmembrane potential) and fix the wave front at time 
t as the hypersurface 

= {x I u{x^t) = Uc and dtu{x^t) > 0}. (7) 

The evolution of the wave front Tip can be described 
by coordinate functions x^ = where a G 

{1,2,...,A^} and A G {1,2, ...,A^ — 1}. Henceforth we 
will use the summation convention for repeated indices. 
Changing from Cartesian coordinates x^ to curvilinear 
coordinates induces a metric 

hAB = dAX°-gabdBx’’ ( 8 ) 



FIG. 1. (Color online) Geometry considered for the varia¬ 
tional principle. The wave front flp (red) with unit normal 
vector ft has passed in all points of the volume V. The in¬ 
terface of V with the domain boundary is a surface Q d with 
outward normal vector ft. Only at the edge C of the wave 
front, the vector field ft is discontinuous; there we define N 
as the outward normal vector to the medium. 


on the wave front with determinant h and matrix inverse 
. In isotropic media, gab is the identity matrix, while 
in systems with anisotropic diffusion it can be taken equal 
to the inverse diffusion tensor to capture spatially varying 
anisotropy [iHiiin]. The coordinate change allows writing 
the linear velocity-curvature relation ^ as: 


= con^ (9) 


with unit normal to the front 
AB...F 

riq = - ^^dAX^dBX\..dFX^ ( 10 ) 

where the e‘ are the fully antisymmetric Levi- 
Civita symbols. Furthermore, in Eq. ([^, = 

{Vh)~'^dA{Vhh^^dB) is n N — 1 dimensional covariant 
Laplacian with respect to the tangential coordinates 
To prove our theorem Q, (|^, we have to show that 


5X 


= n, 



( 11 ) 


under appropriate boundary conditions. When the wave 
front touches the domain boundary, the surface bordering 
the volume where the wave has passed can be divided into 
the wave front surface Clp and the surface Qjj where the 
occupied volume touches the domain boundary. This sit¬ 
uation together with the outer normal vectors is sketched 
in Eig. The N — 2 dimensional border between Qp and 
ftp is C, the wave front boundary. In this work, we will 
not consider broken wave fronts (scroll waves), and as¬ 
sume that wave fronts only end on the domain boundary. 
In this case, the volume where the wave has passed has 
a border Q = ftp U ftp U C that is a closed surface. We 
will assume that the wave front intersects the boundary 
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orthogonally, such that we can impose: 

Vx G 11 n, 

*71 = 0, (12) 

Vx e C : SX is continuous across C. 


The first condition is the transverse gauge, in which we 
convene without loss of generality that points of the wave 
front move orthogonally to it. The second condition en¬ 
sures that the wave front does not propagate across the 
domain boundary, but allows reparameterization of the 
interface. The third condition is compatible with the first 
two only when the wave front ends orthogonally on the 
domain boundary. This situation automatically follows 
if the physical boundary is impermeable to the reactants 
of the underlying system and if the spatial variations in 
the domain boundary are of length scales bigger than the 
wave front thickness. 

In a first step, we construct the volume occupied by the 
wave front. Take the origin of the Cartesian coordinate 
system in the region where the front has passed. The 
volume of the cone with the apex at the origin and an 
elementary surface area of the wave front as the base is 
dV = {1/N)X • fidS with dS = whence 

r AB...F 

V = dS^^€ab..gdAX'^dBX\..dFXfXS. (13) 


We have used the symbol f to emphasize that the in¬ 
tegration runs over the closed surface ft. Under small 
variations 8X{(j^) in the direction of the wave front nor¬ 
mal only, the resulting change of volume is the functional 
differential 

AB..F r 

SV = - f d^-^a[dASX^....dFXJ^6X^ 

• J o 

+ .... + OaX^-.... dpdxf X^ + dAX’^...dFXf 5X^ 

^AB..F^ r 

= - f F~A[-5X‘^....dFXf5dAXs 

- .... - dAX°'....5X^dFXs + dAX‘^...dFXf5XS] 

^AB..F^ p 


= (h dS n-SX. 


(14) 


In the first equality, we performed integration by parts 
and the boundary terms always vanish due to integra¬ 
tion over a closed surface. For the second equality, we 
renamed a GG ^ in the first term, followed by a transposi¬ 
tion in the Levi-Civita symbol, such that the alphabetic 
ordering of indices is restored. Doing this for the (N-1) 
first terms yields the desired result. The last line of Eq. 
( p^ implies that SV/SX = n, proving the first statement 
in Eq. (11). 

Secondly, we vary the wave front surface 


5 ' = 



(15) 


delivering 


SS = 



= - [ d?aVhh^^5hAB 
2Vh 2 J 

d^-^aVhh^^dAX ■ deSX 
d^-AdB{Vhh^^dAX) • 5X + B.T. 
dS V^X ■ 5X. 


(16) 


Therefore, SS/5X = TAX = Kn if the boundary term 
(B.T.) vanishes. Indeed, using the divergence theorem 
with respect to coordinates we find 


B.T. = 


[ d^-AdB{Vhh^^dAX-SX) (17) 

J dp 

f d5div(gradX • (5X) = (f di N ■ 5X = 0. 

J dp J c 


Here, N is the domain boundary at the wave front edge C. 
Then, by Eq. (12), the boundary term vanishes. Hence 
we have demonstrated that relation can be written 
as a gradient system. We conclude that in its motion, 
a non-linear wave front strives to enlarge its occupied 
volume, while minimizing its surface area. 

The corollary © then states that 


‘jS-coV = -clS-'y‘ 


j K^dS+2jco J KdS < 0 . (18) 


which is the analogue of the law for vortex filaments by 
Biktashev et al. [20]. In particular, for dissipative struc¬ 
tures (i.e. when cq = 0 ) which are studied in the context 
of pattern formation, one has that 


S = -l j K^dS. 


(19) 


Thus, for positive surface tension 7 , dissipative structures 
evolve towards a minimal surface. 

Discussion. We have shown that the geometric law 
of motion (§ elegantly follows from an action principle 
that involves only the simplest geometric invariants of 
a closed surface. In physics, variational principles have 
allowed to regard the equations of motion in a more 
abstract manner, and enabled powerful generalizations. 
Eamous examples include the principles of Eermat and 
D’Alembert, the Hamiltonian and Lagrangian approach 
to classical mechanics, and their extensions to quantum 
theory. Other applications of variational principles are 
the approximation of eigenstates using the Ritz method 
and imposing constraints in theoretical and numerical 
problems of fluid dynamics and continuum mechanics. 

Eor non-linear vortices, an extremum principle was 
found for stationary vortex filaments in EUHl]. How¬ 
ever, a front with non-vanishing plane-wave speed cq will 
never reach a final state unless the wave front itself van¬ 
ishes. Therefore, it should be of no surprise that the 
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variational principle which we identified here is of the 
gradient rather than the extremum type. 

At this point, we mention that the terms in the gra¬ 
dient potentials F which we found are highly similar to 
geometric actions in string theory. Namely, the surface 
term (15) is identical to the Nambu-Goto string action 
term, while the volume term (13) has the form of a cou¬ 
pling with an antisymmetric (Kalb-Ramond) field. For 
propagating dissipation waves, this field happens to be 
the volume form of the wave front surface [23]. 


One may wonder whether a similar variational prin¬ 
ciple will hold for vortex filaments in extended media. 
Such filaments can be formed after break-up of a travel¬ 
ing wave front due to the velocity-curvature relation <§• 
A scroll wave filament with arc length parameter s is 
known to obey the geometric law of motion [TSl EO] [24] : 


X(s, t) = -iidlx + -f2dsX X dlx. 


( 20 ) 


If this equation of motion can be written as a gradi¬ 
ent system, it needs to satisfy the compatibility con¬ 
dition Q- However, for the law ( [^ , one computes 


jdX'^ = — 2 j 2 dsX^Cijm^ which does not satisfy the 
criterion unless 72 = 0. This condition is fulfilled in 
chemical systems in which all reactants diffuse at the 
same rate. Thus, the purely dissipative filament dynam¬ 
ics in such systems with no-flux boundary conditions can 
be written as a gradient system with potential 

F = 7 iL, (21) 

where L is the total filament length. The corollary (© 
recovers the known property that filament length changes 
monotonically in time [ 20 ] . 

In summary, we have managed to formulate a simple 
variational principle for the emergent dynamics of propa¬ 
gating non-linear waves in a N-dimensional medium with 
obstacles. Our result is also a contribution to complexity 
theory, since we have been able to identify general laws 
for a wide class of natural systems, without knowing the 
details of the underlying microscopic processes, and even 
without knowing the type of partial differential equations 
that are needed to describe it. 
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